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ABSTRACT 

In this chapter we choose standard airfoil NACA 0015. Which is symmetrical airfoil with a 15% 
thickness to chord ratio was analyzed on ANSYS FLUENT to determine the coefficient of lift, 
coefficient of drag and graph of coefficient of lift v.s\ coefficient of drag. The 2-dimensional cross 
sectional view was considered. The wind velocity was taken as 17m/s which are corresponding to 
232,940 Reynolds number. The airfoil, with an 8 in chord, was analyzed at 0, 5, 10 and 15 degrees. 
Parameters viz. Coefficient of lift (Cl), Coefficient of drag (Cd) and Cl/Cd are calculated and are 
plotted against different angle of attack. 
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1. INTRODUCTION 

NACA AIRFOIL is developed by the NACA (National Advisory Committee for Aeronautics). The shape of 
the AIRFOIL is described by using a series of digits following the word "NACA”. The first digit involves the 
maximum camber as percent of the chord. The second digits involve location of that maximum camber 
measured from leading edge in percent of the chord and last two digits involves maximum thickness o the 
AIRFOIL in percent of chord length. 

The 2-dimensional cross sectional view is considered. The wind velocity is taken as 17m/s which is 
corresponding to 2,32,940 Reynolds number . The angle of attack is varied. 

Parameters Studied - Coefficient of lift (Cl), Coefficient of drag (Cd) are calculated and are plotted against 
different angle ofattack. 

Software - “ANSYS FLUENT” version 13. 1 Jtl 
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NACA airfoils are developed by National Advisory Committee for Aeronautics (NACA). 

Angle of attack 

If you stretch your arm out through the window of car which is moving with some speed, it feels your arm 
pushed backward. When you hold your arm straight with your hand parallel to the road, and change the angle 
slightly, suddenly feel that it is drown upwards. Hand and arm work like the wing of an airplane and with the 
right angle of attack The Computational Fluid Dynamics (CFD) software used is you can feel a strong lift 

[ 3 ] 

force. 

Chord length The distance between leading edge and tailing edge. 

Nomenclature of an Airfoil 



T21 

Figure 1 Basic Nomenclature of an airfoil 


2. METHODOLOGY 

• Pre-analysis & start-up 

• Geometry 

• Mesh 

• Physics setup 

• Numerical solution and results 

2.1. Pre- Analysis & Start-Up 

One of the important things one should think about before one set up the simulation is planning the boundary 
conditions of the given set up. C-Mesh is One of the popular meshes for simulating an airfoil in a stream, and 
that is what we will be using. At the inlet of the system, we have defined the velocity as entering at a 0-degree 
angle of attack and at a total magnitude of 1 . We defined the gauge pressure at the inlet is 0. As for the outlet, 
the only thing we have assumed is that the gauge pressure is 0. As for the airfoil itself, we have treated it like a 
wall. 
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2.2. Airfoil Geometry 



Figure 2 Flow Domain Generation 


Figure 3 Fine meshed flow domain 


[4] 


2.3. Meshing of Flow Domain 

First, we have applied a mapped face meshing control to the geometry. Next, we applied edge sizing [2] 
controls to all of the edges of the mesh. In order to analyze fluid flow, flow domains is split into smaller sub 
domains. The governing equations are applied and solved inside each of these sub domains. The meshing 
around the aerofoil is shown in below figure in which meshing accuracy is increases as we go towards the 
aerofoil. 


Object Name 

Surface Body 

State 

Meshed 

Graphics Properties 

Visible 

Yes 

Transparency 

1 

Definition 

Suppressed 

No 
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Default Coordinate System 

Thickness 
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Material 

Fluid/Solid Defined By Geometry (Solid) 
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3. RESULTS AND DISCUSSIONS 

In this study, numerical analyses were performed. The analysis was performed at 17 m/s wind velocity. 
Coefficient of Lift and drag of NACA 0015 airfoil at different angle of attack between 0° and 15° were 
calculated. The lift and drag coefficients are obtained as numerical with FLUENT programs for the same 
conditions. In numerical analysis C mesh was used. The top bottom and left boundaries are placed at a 
distance of 10 chords from airfoil. Whereas the right boundary was placed at 20 chords. A mesh independent 
study was performed to verify that the solution would not change subsequent additional refinements and 
1,35,000 grids number suitable for our model. Airfoils have various shape and sizes. The non- dimensional 

coefficients for two dimensions’ solution is given as below: ^ 

Lift coefficient [6] : 

Cl = 2L/gv 2 c 
Drag coefficient [6] : 

Cd =2D/gav 2 c 

Where, L and D are lift, and drag force, Cl and Cd are lift and drag coefficient of airfoil respectively, c is 
airfoil cord length, V is velocity of wind, p is density of air. 

3.1. Results 

The lift coefficient varies with angle of attack. Increasing angle of attack is associated with increasing lift 
coefficient. After maximum lift coefficient, lift coefficient decreases as the angle of attack. A symmetrical 
wing has zero lift at 0 degrees angle of attack. 

Identically the value of drag coefficient is zero at the zero angle of attack and it increase slowly till the stall 
condition and at the time of stall as well as after stall it increase readily. 

The critical or stalling angle of attack is typically around 15° for many airfoils. 

Coefficient of drag and lift with iterations at different angle of attack 
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Figure 4 Plot between coefficient of drag and iterations at an angle of attack 0 degree 
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0 200 400 600 600 1COO 1200 1400 1600 1000 2000 

Iterations 

Figure 5 Plot between coefficient of lift and iterations at an angle of attack 0 degree 



Figure 6 Plot between coefficient of lift and iterations at an angle of attack 5 degree 



Figure 7 Plot between coefficient of lift and iterations at an angle of attack 10 degree 
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Figure 8 Plot between coefficient of drag and iterations at an angle of attack 15 degree 

3.2. Velocity Counters 
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Figure 9 Velocity counter at an angle of attack 0 degree 


■ 


2.580+01 
2.450+01 
2 , 330+01 
2 . 210+01 
2,Q9*+01 
1 .960+01 
1.840+01 

1 720+01 
1 590+01 
1 .470+01 
1 35et01 
1 23.+01 

1 -IO 0 +OI 
9.79e+O0 
8.56e+O0 
7.340+00 
6 . 110+00 
4.88e+O0 
3.650+00 

2 420+00 
1 190+00 





Figure 10 Velocity counter at an angle of attack 5 degree 
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Figure 11 Velocity counter at an angle of attack 10 degree 
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Figure 12 Velocity counter at an angle of attack 15 degree 


3.3. Pressure Counter 

At zero lift the pressure distributions over the upper and lower surfaces is identical. A contour of static 
pressure shows that static pressure increases at the lower surface of the aerofoil with increasing angle of 
attack. 

At no incidence, In case of symmetrical airfoil, the distribution of velocity and thus the pressures along 
both surfaces would have been exactly the same, canceling each other & total lift force is zero. 
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Figure 13 Pressure counter at an angle of attack 0 degree 
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Figurel4 Pressure counter at an angle of attack 5 degree 
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Figure 15 Pressure counter at an angle of attack 10 degree 



Figure 16 Pressure counter at an angle of attack 15 degree 
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Figure 17 Cd vs angle of attack 



Figure 18 Cd vs Cl 
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Figure 19 Cl vs angle of attack 
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Table 1 Comparison of Theoretical and calculated values of properties [1][5] 





Calculated 
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- 

16 

- 

- 
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4. CONCLUSIONS 

In this lift and drag performances of NACA 0015 airfoil were performed. A ANSYS FLUENT program was 
used to numerical calculations. Numerical and experimental results were compared. The inference from 
calculated results is as follows: 

• Drag and lift coefficients increased with increasing angle of attack. 

• Stall was started with 16° attack angle. Lift coefficient increased whereas; drag coefficient increased. 

• The optimum lift coefficient value was computed atl6° 
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